Sine-square deformation of free fermion systems in one and higher dimensions 
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We study free fermion systems with the sine-square deformation (SSD), in which the energy scale 
of local Hamiltonians is modified according to the scaling function f{x) = sin^ [x'(^ ~ ^h^re 
X is the position of the local Hamiltonian and L is the length of the system in the x direction. 
It has been revealed that when applied to one-dimensional critical systems the SSD realizes the 
translationally-invariant ground state which is the same as that of the uniform periodic system. In 
this paper, we propose a simple theory to explain how the SSD maintains the translational invariance 
in the ground-state wave function. In particular, for a certain one-dimensional system with SSD, it 
is shown that the ground state is exactly identical with the Fermi sea of the uniform periodic chain. 
We also apply the SSD to two-dimensional systems and show that the SSD is able to suppress the 
boundary modulations from the open edges extremely well, demonstrating that the SSD works in 
any dimensions and in any directions. 

PACS numbers: 71.10.Fd, 71.10.-w, 05.30.Fk 



I. INTRODUCTION 

The boundary condition has a crucial influence on 
properties of quantum systems. For example, the peri- 
odic boundary condition (PBC) realizes a translationally- 
invariant system with the geometry of the ring or torus, 
while the open boundary condition (OBC) results in 
the chain or cylinder with open boundaries with the 
boundary-induced modulation such as the Friedel os- 
cillations. The thermodynamic properties, i.e., the ex- 
pectation values of observables in the bulk of infinite- 
size systems, are believed to be the same irrespective of 
the boundary conditions employed. However, the way 
that the expectation values approach the thermodynamic 
limit is usually distinct for OBC and PBC, which requires 
a different treatment of finite-size scaling. Furthermore, 
at the level of the wave function, the difference between 
the critical systems under OBC and PBC remains siz- 
able even at the thermodynamic limit. That is clearly 
demonstrated by the different scaling behavior of the en- 
tanglement entropy in one-dimensional (Id) critical sys- 
tems with and without open boundaries ii 

Attempts to control the boundary and finite-size effects 
have been made in the studies of lattice fermion/ spin 
models. One example is the so-called smooth bound- 
ary condition, which succeeded pretty well in suppress- 
ing the boundary effects by smoothly decreasing the 
energy scales near the edges of the system.—— Similar 
but modified boundary conditions were applied to trans- 
port problems such as the conductance of Id interacting 
systemsi^ Another example is the hyperbolic deforma- 
tion in which the ground state expectation values of local 
observables are nearly uniform in the bulk whereas the 
interaction strength is enhanced rather than suppressed 
at the boundaries)^"— 

Recently, a more efficient scheme to suppress the 
boundary effects, which is called the sine-square defor- 



mation (SSD), has been proposed for Id critical systems. 
In the system with SSD, the energy scale of the local 
Hamiltonian is rcscalcd according to the scaling function. 
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where L is the length of the system and x is the center po- 
sition of the local Hamiltonian. Note that since f{x) = 
for X = 1/2 (mod L), the SSD disconnects the link be- 
tween the sites at a; = 1 and L, resulting in the open 
edges. The SSD was first applied to the Id free-fermion 
system.— It was found there that in the system with SSD 
the local observables such as the bond strength are trans- 
lationally invariant, suggesting that the SSD was able 
to remove the open-boundary effects almost completely. 
Furthermore, the effects of SSD in several quantum spin 
models at criticality were examined numerically.— In ad- 
dition to the translational invariance of the observables 
such as spin correlations, it was further revealed that the 
wave-function overlap between the ground state of the 
open system with SSD and that of the uniform system 
with PBC is very close to (almost exactly) unity. This 
means that the SSD docs not disturb the periodic ground 
state at the level of the wave function and keeps the state 
as the true ground state even after the model Hamilto- 
nian loses the translational symmetry by the energy de- 
formation. This striking feature of the SSD was proven 
rigorously for the spin- 1/2 XY chain, which is equivalent 
to the free fermions in one dimensionJ^ It has also been 
shown that the SSD works well for other strongly cor- 
related systems including the Id Hubbard modelii and 
Kondo-lattice modelJ^ 

While the efficiency of the SSD has been confirmed for 
several Id systems, it remains a puzzle why the ground 
state of the system with SSD is (almost) identical to that 
of the periodic system. Since the system with SSD is no 
longer translational invariant, its one-particle eigenstatcs 



are distinct from the plane waves. However, when (and 
only when) the fermions are filled up to the Fermi level, 
the many-particle ground states of the system with SSD 
and the uniform periodic system become equivalent. It 
is desirable to clarify the mechanism of the SSD to leave 
the ground state under PBC unchanged. 

In this paper, we propose a simple theory which ex- 
plains how the SSD realizes the energy-scale deformation 
with keeping the ground state of the original free-fermion 
model with PBC as an eigenstate. The key idea is that 
the amplitudes of active processes in the Hamiltonian 
with SSD vanish, therefore the Hamiltonian does not dis- 
turb the Fermi sea of the original model with PBC. For 
a certain class of Id systems, it is also proved that the 
Fermi sea of the original uniform system is not only an 
eigenstate but also the unique ground state of the Hamil- 
tonian with SSD. For two-dimensional (2d) systems, the 
ground state of the system with SSD is found to be very 
close to the Fermi sea of the original periodic system ac- 
companied by edge states localized around open edges. 
Wc emphasize that the SSD works well to suppress the 
boundary effect induced by the OBC in any dimensions 
and in any directions. The SSD thereby reduces the 
number of directions under PBC with maintaining the 
ground-state properties in the bulk. For example, apply- 
ing the SSD to a 2d model changes its geometry from a 
torus to a cylinder, then to a rectangle. 

The paper is organized as follows. In Sec. [Til we discuss 
the Hamiltonian with SSD in the momentum space and 
show that the Fermi sea of the uniform system remains 
an (approximate) eigenstate through the energy defor- 
mation. The results for the Id and 2d systems are pre- 
sented in Sec. Ill A] and III Bl respectively. In Sec. IIIII we 
study numerically the effects of the SSD on the energy- 
level structure as well as the ground state properties of 
2d systems. By investigating several quantities such as 
the one-particle eigcn-cnergies, eigcn-wavefunctions and 
the ground-state density profiles, we examine how the 
ground state is modified by the SSD. We conclude with 
some remarks in Sec. IIVI 



II. SSD IN MOMENTUM SPACE 

A. One-dimensional systems 

In this section, we develop a theory of the mechanism 
of the SSD applied to free spinless fermions on a lattice. 
We first focus on the simplest model, i.e., the Id tight- 
binding model with only nearest-neighbor hopping. For 
the model, it was rigorously proved that the model with 
SSD shares the same ground state with the model under 
PBC^ 

The model Hamiltonian of the tight-binding chain is 
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FIG. 1: Dispersion curve and Fermi sea of the Id tight- 
binding model ([2]) with t > and ^ < 0. Dark (blue) and 
light (green) arrows represent respectively the processes in 
the chiral Hamiltonian H^^^ and "H^"' allowed by the Pauli 
principle. 



given by 

H = -t^[c\x)c{x + 1) + R.C.] - fi^c'ix)c{x), 

X X 

(2) 

where c^{x) and c{x) are respectively the creation and an- 
nihilation operators of fermion at the site x. We consider 
the L-site system and impose the PBC. The Hamiltonian 
in momentum space is obtained as 

?i = ^e(/c)4cfe (3) 

k 

with the dispersion relation 

e{k) ^ ~2tcosk ~ fi. (4) 

A key ingredient of the construction of SSD is to in- 
troduce the chiral deformation, which leads to the chiral 
Hamiltonian of the form, 

H^^'> = -t^e±*^" [ct(x)c(a; + l) + H.c.] 

X 

-t,Y,e^^'<--'i^c\x)cix), (5) 

X 

where S ~ 2tt /L. The chiral Hamiltonian 'H^^'^ is nothing 
but the Fourier component of the local Hamiltonians and 
is apparently non-Hermitiauji^ In the basis where the 
original Hamiltonian % is diagonal, H*^^^ is represented 
as 

H(±) = e^'hik T S/2)clck^s. (6) 
fe 

The point to be noted here is that the chiral Hamil- 
tonian includes only the terms with momentum trans- 
fer Afc = ±(5, which connects the nearest-neighbor sites 
in k space, and their amplitude is proportional to the 
dispersion function Eq. Q of the original uniform sys- 
tem. Hence, when applied to the Fermi sea of the original 
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model, most of the transfer processes are prohibited by 
the Pauh principle; Only the processes across the Fermi 
point k = ifcp can be active. (The schematic picture 
of the active processes is shown in Fig. [TJ) However, for 
such processes, the amplitudes are, at most, of the or- 
der of 0{1/L). Indeed, one can make these amplitudes 
be exactly zero by fine-tuning of the chemical potential 
fi so that e(fc =F i5/2) = for the active processes. This 
means that the Fermi sea of the original uniform chain is 
an exact eigenstate of the chiral Hamiltonian 'H'^^'> and is 
not disturbed by the chiral deformation at all. The same 
result holds for any Id systems having a single Fermi 
momentum fcp and symmetric dispersion e(fc) = e(— fc). 

Now, we are ready to consider the system with SSD. 
The Hamiltonian with SSD is constructed from the orig- 
inal and chiral Hamiltonians as 

= -t^/(a; + l/2) [ct(a;)c(x + l)+H.c.] 

X 

-Ai^/(x)ct(x)c(x), (7) 



where the scaling function f{x) is given by Eq. ([T]). Here, 
the chemical potential fi has been chosen in such a way 
that the Fermi sea of Ji for a fixed number of particles is 
annihilated by H*-^-*. Since f{x) = for x = 1/2 (mod 
L), the links between the sites with x = 1 and L are dis- 
connected by the SSD, resulting in the open boundaries 
at the edges x = 1 and L. From the argument for the 
chiral Hamiltonians above, it follows that the Fermi sea 
of the original periodic model is an exact eigenstate of 
the Hamiltonian with SSD. 

So far, we have seen that the Fermi sea of the uni- 
form periodic system is an exact eigenstate of the system 
with SSD. For the model with the only nearest neigh- 
bor hopping, one can further show that the Fermi sea is 
indeed the exact ground state of the SSD Hamiltonian. 
The proof is essentially the same as the one presented in 
Rcf. ilOi and we briefly review it here. We first transform 
both the original and SSD Hamiltonians into the XY spin 
chains using the Jordan- Wigner transformation. Then, 
for fixed magnetization (for fixed number of fermions), 
it can be shown that (i) all of the off-diagonal elements 
of the Hamiltonian are non-positive, and (ii) the Hamil- 
tonian satisfies the connectivity condition, for both the 
PBC and SSD cases. Therefore, the Perron-Frobenius 
theorem applies and guarantees that the ground state is 
unique in each case. Since the Perron-Frobenius theorem 
tells us that the ground state is nodeless, i.e., can have 
only positive components in an appropriate basis, the 
Fermi sea of the PBC Hamiltonian is the ground state 
of the SSD Hamiltonian if it has the Fermi sea as an 
eigenstate. 



B. Two and higher dimensions 

In the preceding section, we have shown for the Id 
free fermion system that the chiral and SSD Hamiltoni- 
ans have the Fermi sea of the uniform periodic system 
as an exact eigenstate. In this section, we extend the 
theory to 2d or higher-dimensional systems. We focus 
on the single-band model, having single site in each unit 
cell. (The multi-band case will be discussed in the end 
of this section.) We employ the 2d square lattice as a 
typical example, but the extension to the other lattices 
and dimensions is straightforward. 

Let us consider the tight-binding model in the square 
lattice, 

n = -tJ2[c\x:y)c{x + l,y) + }l.c.] 

r 

-t_L 2/ + 1) + H.C.] 

r 

-tiJ2[cHx,y)c{x + l,y + l)+R.c.] 

r 

-t2 J2 y)<^ + l,y-l)+ H.c] 

r 

-fJ'Yc^{x,y)c{x,y). (8) 

r 

We set the lengths of the system along x and y directions 
Lx and Ly, respectively, and impose the PBC in both 
directions. In addition to the nearest-neighbor hoppings 
t and tj_ in x and y directions, we have included the 
diagonal hoppings ti and t2 for generality. The schematic 
picture of the model is shown in Fig. [5] (a). In momentum 
space, the Hamiltonian is diagonalized as 

n = Y.e{k)cick (9) 
fc 

with the dispersion relation 

e(fc) = — 2t cos /cj; — 2tj^ cos fcj, 

— 2ti cos(fc2, + ky) — 2t2 cos{kx — ky) — fi, 

(10) 

where fc = (k^, ky). 

Similarly to the Id case, one can introduce the chiral 
deformation to the system. The resulting Hamiltonian is 

= -t^e±*[*=^+*''(''-^)l [ct(a;,y)c(.T + l,y) + H.c.] 

r 

-t±Yl e^'^^^-^^'-i^+^yy^ [c^ix, y)c{x, y + 1) + H.c] 

r 

-tiJ2 e^'^^-'^+^^y^ [c^x, y)c{x + l,y + l)+ H.c] 

r 

- ^2 XI e^'^^^^'+^y^y-^'^'^ [ct(x, y)c{x + 1, y - 1) + H.c] 

r 

- X e±''[^-("-^)+*«(^-^)lct(x, 2/)c(x, y), (11) 
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FIG. 2: (a) Schematic picture of the square-lattice model 
(b) Fermi sea of the model (O with t = t_L, ti = t2 = 0, 
and fi < 0. Dotted line shows the Fermi surface. Dark (blue) 
and light (green) arrows represent respectively the processes 
allowed by the Pauli principle in the chiral Hamiltonian 
and n'f'' with <5 = {2n/L^,0). 



where we consider the cases of <5 = {Sx,Sy) = 
{27r/Lx,0),{0,2n/Ly), or {2'K/L^,±2n/Ly). From the 
usual Fourier transformation, the chiral Hamiltonian in 
k space is obtained as 

^i^' =E^^^'^'^^'''^<fcT^/2)cU^6. (12) 

k 

Hence, we find again that the chiral Hamiltonian contains 
only the processes connecting the sites k and fc ± ^ with 
an amplitude proportional to the dispersion function Eq. 
()10p . Among them, the processes being active are only 
those across the Fermi surface, whose amplitudes are of 
the order of ©(l/L^) or 0{1/Ly). Therefore, the chiral 
Hamiltonian does not disturb the Fermi sea very much; 
The effect of the mixing is of the order of 0{1/Lx,y) and 
vanishes in the limit L^^y — s- oo. 

Here, the difference from the Id case should be noted. 
For Id systems, one can make all the amplitudes of the 
disturbing processes be exactly zero by tuning the chemi- 
cal potential fi as discussed above, and then the Fermi sea 
can be an exact eigenstate of the chiral Hamiltonian even 
at finite or Ly. For systems in higher dimensions, on 
the other hand, it is impossible to achieve e{k^S/2) = 
for all the processes on the Fermi surface by tuning a sin- 
gle parameter fj,, except for some cases in which the Fermi 
surface takes a special shape. Therefore, the Fermi sea is 
an approximate eigenstate of the chiral Hamiltonian and 
the correspondence becomes exact only asymptotically 
in the limit L^.y — > oo. Wc note that one can achieve 
the exact correspondence at finite and/or Ly by in- 
troducing momentum-dependent chemical potential. For 
example, when the chiral deformation in the x direction 
with S = {2tt/ Lx,0) is applied, one may introduce ky- 
dependent chemical potential fJ-{ky), which leads to long- 
range hoppings along the y direction in real space, and 
fine-tune them to realize the zero amplitude for all active 
processes in the chiral Hamiltonian. 

Another notice for higher dimensions concerns the or- 
der of taking the infinite-size limits in the x and y 



directions. For example, let us consider the case of 
S = {2Tr/Lx,0)- In this case, the amplitudes of the ac- 
tive processes are of the order of 0(1/Lx) and vanish at 
Lx — > oo. However, since the number of these processes 
is proportional to Ly, their effect may be amplified as 
Ly increases. In the argument of this section, we take 
the limit L^ ^ oo first then consider large enough Ly. 
Taking the limit L^ = Ly ^ oo simultaneously is rather 
subtle, however, the numerical results presented in Sec. 
Hill suggest that the SSD works well also for the case of 

Lx = Ly. 

Let us move to the system with SSD. By taking an 
appropriate linear combination of the original and chiral 
Hamiltonians, one can construct various kinds of the SSD 
Hamiltonian, which has the form, 

= -t^J'ix + ^^y) [c\x,y)c{x + l,y) + ll.c.] 

r 

-t±Y^ T{x, y + i) [c^(x, y)c{x, y + l)+ H.c] 

r 

- E ^(a; + i, y + i) [c\x, y)c{x -M, y + 1) + H.c] 

r 

-t2Y,H^ + \^y-\) Wi^^ y)<^' + 1, y - 1) + H.c] 

r 

~ ^i^J'{x,y)c){x,y)c{x,y). (13) 

r 

The scaling function J-{x, y) can take various forms; for 
example, F{x,y) = fx{x) for the SSD only in the x di- 
rection and J-{x,y) = fx(x)fy{y) for the SSD in both 
the X and y directions, where fx(x) [fy{y)\ is given by 
Eq. dH) with L Lx {Ly). It is clear that these SSD 
Hamiltonians have the Fermi sea of the original model as 
an approximate eigenstate. 

We have shown that the chiral and SSD Hamiltonians 
have the Fermi sea of the uniform periodic system with 
small disturbance of 0{1/Lx^y) as an eigenstate. How- 
ever, the argument docs not tell us whether the Fermi 
sea is the ground state of the system with SSD. As men- 
tioned in the preceding section, for the Id system, the 
Perron- Frobenius theorem is applicable and enables us to 
prove that the Fermi sea is the ground state of the SSD 
Hamiltonian. In contrast, in higher dimensions, it is in 
general impossible to find a basis in which all of the off- 
diagonal elements of the Hamiltonian are non-positive, 
due to fermion signs. Therefore, the ground state of the 
SSD Hamiltonian is not necessarily the Fermi sea of the 
periodic system, and one must examine how the ground 
state evolves by the SSD using some other methods. We 
will show in Sec. IIIII numerically that in two dimension 
the ground state of the system with SSD is not exactly 
the Fermi sea but its slight modification accompanied by 
edge modes locahzed at boundaries. 

The result above is valid regardless of the boundary 
conditions in the direction to which the SSD is not ap- 
plied. For example, when applying the SSD in the x 
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FIG. 3: (a) Energy levels of the one-particle eigenstates of 
Hia) for t = t± = l,ti = t2 = fi — 0, and — Ly = 7 
and (b) its enlarged view around a = 1. In (b), thin-dashed 
and bold lines represent single and doubly-degenerate levels, 
respectively, (c) Fidelity between the Fermi sea at a = and 
the many-particle ground state of 'H(a) at a > with N = 25 
particles, (d) Normalized density profiles LxLy\ipi{x,y)\^ of 
25, 26, and 27th one-particle eigenstates at a = 1 as a function 
of X. Note that the profiles are independent of y because of 
the PBC imposed in the y direction. 



direction, one can develop the same logic for the system 
with OBC in the y direction as long as the magnitude 
of the momentum of the direction is a good quantum 
number. The SSD thus changes the boundary condition 
in a direction from PBC to OBC with the energy de- 
formation, leaving the boundary conditions in the other 
directions unchanged. 

An interesting question is whether one can extend the 
above argument to the multi-band cases where there are 
more than one site per unit cell. Unfortunately, it turns 
out that the chiral deformation gives rise to interband 
hopping terms in the dual lattice (fc-space) in these cases. 
Since the Pauli principle does not prohibit the interband 
hopping processes, the Fermi sea is not expected to be an 
approximate eigenstate of the SSD Hamiltonian except 
for some special cases. Therefore, we can conclude that 
the validity of the SSD is basically limited to single-band 
systems. 



III. NUMERICAL RESULTS 
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FIG. 4: Numerical results for the square-lattice model with 
diagonal hoppings, t — t± = l,ti = t2 = 0.5, n = 0.3, and 
Lx — Ly = 60; (a) Density of state, where we use the Gaus- 
sian broadening with a — 0.5. (b) Fidelity between the Fermi 
sea at a = and the many-particle ground state at a > 
with fixed particle number, (c) Normalized density profiles 
LxLy\(pi{x,y)\^ at a = 1. The results for eleven states with 
the smallest \Ei \ are plotted, (d) as a function of \Ei\. The 
data are normalized by the value for the uniform periodic 
system, A?(a = 0) = {L^ - 1)(L^ - 2)/12. The horizontal 
dotted line is a guide for the eye. 



In this section, we discuss how the ground state 
changes when we apply the SSD to the 2d system under 
PBC. For the purpose, we introduce the model connect- 
ing the uniform system and the system with SSD, 



n{a)^{l-a)n + aHssD, (14) 
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where we consider the case that the SSD is apphed in x 
direction. Therefore, T-L{a) is given by the r.h.s. of Eq. 
(fTS]) with replacing the scaling function J-{x, y) by 

/(x;a) = (l-a) + a/,(x). (15) 

Note that a ~ corresponds to the uniform system with 
PBC while a = 1 gives the system with SSD. 

First, we consider the simple square lattice with = 
Ly ^ L, t = t± = 1 and ti = t2 = H = 0, where 
the Fermi surface has the diamond shape for infinite L. 
For finite and odd L, the amplitude e{kx T ^jky) in 
the chiral Hamiltonian, Ea. ((T^ . becomes exactly zero 
for the processes across the Fermi surface, therefore the 
Fermi sea of the uniform model with a = is an exact 
eigenstate for arbitrary a. For even L, the Fermi sea 
at a = becomes an exact eigenstate for arbitrary a 
when the boundary condition in the y direction is anti- 
periodic. We again emphasize that an exact eigenstate 
does not mean the ground state. In fact, the Fermi sea at 
a = is, in general, not the ground state but an excited 
state for a = 1 due to the level crossing, as shown below. 

Figures [31Ja) and (b) show the energy levels of the one- 
particle eigenstates of ?^(a) for L = 7 with PBC in the 
y direction. With increasing a, many energy levels are 
approaching to i? = 0. At a — 0.982, a state crosses the 
Fermi energy E = 0, leading to the change of the particle 
number of the ground state filling the negative energy. 
Furthermore, at a = Uc = 0.996, the level crossing occurs. 
Therefore, even if the particle number N is fixed at that 
of a = 0, the ground state changes abruptly at a = Oc 
and the Fermi sea at a = is no longer the ground state 
for a > Qc- To elucidate it, we plot in Fig. ^c) the 
fidelity, i.e., wave- function overlap, between the Fermi sea 
at a = and the ground state of 'H(a) for fixed = 25 
particles. The fidelity jumps from 1 to at the level 
crossing point a = ac. We also confirmed numerically 
that the Hamiltonian ?^(a) has the Fermi sea at a = as 
an exact, first-excited eigenstate for a > ac- 

Here, it is non-trivial that the ground state at a = 
defined by the Slater determinant of the negative-energy 
states is the exact eigenstate of TL{a) in spite that the 
one-particle eigenstates are different. While the local 
density of one-particle states at a = 0, which arc a sim- 
ple plane wave, does not depend on the spatial position 
{x,y), the density for a > depends on x due to the 
energy deformation in x direction. In Fig. [3l^d), we show 
the normalized density profiles LxLy\ipi{x,y)\'^ of three 
low-energy one-particle states ipi{x,y) at a = 1, which 
states undergo the level crossing at a = Cc. As shown 
in the figure, those states are localized at the boundaries 
at a; = 1 and L. The result implies that the difference 
between the Fermi sea of the uniform system H = T-L{0) 
and the ground state of Hssd = "^(1) manifests itself in 
the boundary region. 

Next, we show that the above features of SSD arc 
commonly seen in the more generic case t =^ t± = 1, 
ti = ^2 = 0.5, and fi = 0.3, where the Fermi sea of 
7^(0) is not the exact ground state nor even an exact 



eigenstate for a > due to small disturbance of 0{1/Lx) 
by the chiral Hamiltonians. Figure |H[a) shows the den- 
sity of states at a = and a ~ 1 for the systems with 
Lx = Ly = 60. The prominent peak at E = for a = 1 
indicates that many levels approach E = after the SSD, 
as in Fig. ^a). 

Figure IH^b) shows the fidelity between the Fermi sea 
at a = and the ground state at a > with the fixed 
particle number. Due to the small disturbance mentioned 
above, the fidelity is not exactly one even for small a in 
contrast to Fig. ^h). In the usual perturbation theory 
via small a, the fidelity is linear in a and its gradient is 
determined by the expectation value of the perturbative 
Hamiltonian in the unperturbed ground state. Therefore, 
the fact that the fidelity at small a remains almost unity 
means that the expectation value of 'H^^-' is very small. 
The fidelity exhibits a continuous decrease from unity at 
a = to almost zero at a = Oc = 0.9991, followed by a 
small jump to zero at a = Uc [which is almost invisible in 
Fig.HJb)]. We have observed that the finite jump of the 
fidelity, which is due to the level-crossing of the ground 
state and excited state, generally occurs at a = Oc ~ 1. 

We present in Fig. UJ^c) the normalized density profiles 
of the one- particle states with smallest \Ei\ at a = 1, 
which are expected to be the most responsible for the 
change in the ground-state properties. We note that the 
figure plots eleven states locating closest to the Fermi en- 
ergy. It is clear that the profiles of those states are almost 
identical and strongly localized around the open bound- 
aries X = 1 and L^, in marked contrast to the plane waves 
for a = whose profiles are a constant. We note that the 
decay of the profiles from the open edges seems to be al- 
gebraic, while determining the detailed form of the decay 
function as well as the decay exponent requires a more 
systematic analysis, which is left for future studies. 

To further clarify the edge-localized nature of the low- 
energy eigenstates, we introduce a measure of the local- 
ization, 

x,y 

Note that A? is small for localized states while it be- 
comes large for the states extended to the bulk. For the 
uniform system H = 'H(O), Af is easily obtained to be 
Af = {Lx - l){La: - 2)/12 for even L^, since \ipi{x,y)\'^ = 
l/{LxLy) for all eigenstates. For a = 1, we obtain numer- 
ically all of the eigen-energies and eigen-wavefunctions 
{Ei,ipi) and calculate A|. Figure Sljd) indicates that all 
the eigenstates with small \Ei \ have very small A| This 
means that the states around the Fermi level, which get 
closer and converge to = as a increases from to 1, 
are strongly localized around the open edges. The local- 
ization of low-energy states can be naturally understood 
from the fact that the scaling function of the SSD, fx{x), 
is small around the open edges, therefore the excitations 
around the edges can appear with an extremely low en- 
ergy cost. 
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FIG. 5: Density profiles p{x, y) in the ground state of (a) the 
square-lattice model, t = t± — 1, ti = t2 — 0, and jj, = 0.3, (b) 
the square-lattice model with diagonal hoppings, t = t± = 1, 
ti — t2 = 0.5, and fi — 0.3, and (c) the triangular-lattice 
model, t — t± = ti = 1, t2 = 0, and fi — 0. The system size 
is = Ly = 60. Solid and open circles respectively show 
the results for the system with SSD and the uniform system 
with OBC in the x direction. The PBC is imposed for the y 
direction, therefore p{x, y) is independent of y. Insets present 
the enlarged figure of the same data for 15 < x < 46. 



Finally, we investigate how the SSD affects properties 
of the many-body ground state. In Fig. [Sj we show the 
density profiles 

p{x,y) ^ {c){x,y)c{x,y)) = ^ \ip^{x,y)\^ (17) 
in the ground state of the system with SSD in x direc- 



P(x.y) 




P(x.y) 




FIG. 6: Density profiles p{x,y) in the ground state of the 
square-lattice model, t = t± = 1, ti = t2 = 0, p = 0.3, 
and Lx ~ Ly = 60; (a) the result for the OBC imposed in 
both X and y directions, and (b) for the SSD in both the 
directions. The data for 30 < y < 60 are plotted while those 
for 1 < y < 30 are symmetric with respect to y — {Ly + l)/2. 



tion [Eq. with T{x,y) = fx{x)] for several typical 
parameter sets. The density y) in the ground state 
of the uniform system with OBC in x direction is also 
shown for a comparison. Remarkably, the density profile 
p{x^ y) in the systems with SSD is almost a constant in 
the bulk and the Friedel oscillations pronounced in the 
uniform system with OBC are strongly suppressed. The 
steep change in p{x, y) from the constant is found only in 
the vicinity of the open edges. The results indicate that 
the SSD preserves the translationally-invariant nature of 
the Fermi sea of the periodic ground state in the bulk 
and the deviation appears only around the open edges. 

Figure |6] shows the density profiles in the square- 
lattice system with SSD in both x and y directions, 
whose Hamiltonian is given by Eq. (jl3p with J-{x, y) — 
fx{x)fy{y)- The result for the system with OBC in both 
the directions is also shown. Again, we find that p{x, y) 
is nearly a constant in the bulk of the system with SSD 
and the Friedel oscillations, which are conspicuous in the 
system with OBC, are removed almost completely. In 
the system with SSD, the deviation in p{x, y) from the 
constant is sizable only around the open edges as well as 
the corners. We have observed essentially the same re- 
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suits for the square-lattice model with diagonal hoppings 
and triangular-lattice model. The completely flat profiles 
strongly suggests that the ground state of the systems 
with SSD is very close to the translationally-invariant 
ground state of the uniform periodic system superposed 
by the small disturbance of the edge states. The SSD 
thereby realizes the topology change of the system from 
a torus to a cylinder, and then, to a rectangle with keep- 
ing the ground-state properties in the bulk unchanged. 

IV. CONCLUDING REMARKS 

We have studied effects of the SSD on the free fermion 
systems. We have proposed a simple theory of the mech- 
anism of the SSD to realize the energy deformation with 
keeping the ground state of the periodic system un- 
changed. For a certain class of Id systems, it is shown 
that the Fermi sea of the uniform system remains an ex- 
act ground state of the system with SSD. For 2d systems, 
it is numerically found that the SSD realizes almost fiat 
density profiles in the bulk, suggesting that the ground 
state properties in the bulk are translationally invariant 
in spite that the system has open edges. The ground 
state of the two- and higher-dimensional systems with 
SSD can be regarded as a Fermi sea of the uniform peri- 
odic system superposed by edge-localized states. 

As we have shown, in one dimension, the Fermi sea is 
not only an exact eigenstate but also the ground state 
of the Hamiltonian with SSD if the chemical potential 
fi is chosen so that e(fc ^ d/2) = for the processes 
in the SSD Hamiltonian allowed by the Pauli principle. 
One can explicitly confirm this correspondence in a class 
of exactly solvable spin chainsJ^ It is worth mentioning 
that the Hamiltonian with SSD is distinct from the in- 
homogeneous integrable models obtained from the quan- 
tum inverse scattering method and the algebraic Bethc 
ansatz ji^iii From the point of view of conformal field the- 



ory, the chiral Hamiltonian has a definite meaning, which 
accounts for the correspondence between the PEC and 
SSD Hamiltonians in a wide class of Id critical systems. 
A detailed discussion will also be given in a subsequent 
paper 

Our results demonstrate that the SSD works ex- 
tremely well for free fermions in any dimensions. A 
natural question to ask is whether the results gener- 
alize to systems with interactions. In one dimension, 
it has been confirmed numerically for several strongly- 
correlated fermioniiii^ and spin^ systems that the SSD 
works efficiently to realize the translationally-invariant 
ground state, and the results are well founded by a field- 
theoretical analysisJ^ Extending the analysis to multi- 
leg ladder systems^i^ and then, to 2d systems must be 
an intriguing problem. It would also be interesting to 
explore the generalizations of other deformations such as 
the hyperbolic^'"— exponential)^ and sinusoidalii defor- 
mations, to higher-dimensional and/or interacting sys- 
tems. 

An experimental study of the SSD in real systems 
should also be a challenge of great interest. Considering 
the fact that hopping strength and coupling parameters 
can be tuned in systems of ultracold atoms in optical lat- 
tices as demonstrated in Refs.[20l-[2^ they can be promis- 
ing candidates for the realization of the system with SSD 
studied in our work. 
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